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The inspiral of a "small" ~ 1 - 100 Mq) compact body into a "large" (M ~ 10^"^ Mq) black 
hole is a key source of gravitational radiation for the space-based gravitational-wave observatory 
LISA. The waves from such inspirals will probe the extreme strong-field nature of the Kerr metric. 
In this paper, I investigate the properties of a restricted family of such inspirals (the inspiral of 
circular, inclined orbits) with an eye toward understanding observable properties of the gravitational 
waves that they generate. Using results previously presented to calculate the effects of radiation 
reaction, I assemble the inspiral trajectories (assuming that radiation reacts adiabatically, so that 
over short timescales the trajectory is approximately geodesic) and calculate the wave generated as 
the compact body spirals in. I do this analysis for several black hole spins, sampling a range that 
should be indicative of what spins we will encounter in nature. The spin has a very strong impact 
on the waveform. In particular, when the hole rotates very rapidly, tidal coupling between the 
inspiraling body and the event horizon has a very strong influence on the inspiral time scale, which 
in turn has a big impact on the gravitational wave phasing. The gravitational waves themselves are 
very usefully described as "multi- voice chirps": the wave is a sum of "voices", each corresponding 
to a different harmonic of the fundamental orbital frequencies. Each voice has a rather simple phase 
evolution. Searching for extreme mass ratio inspirals voice-by-voice may be more effective than 
searching for the summed waveform all at once. 

PACS numbers: 04.30.Db, 04.30.-w, 04.25.Nx, 95.30.Sf 



I. INTRODUCTION 

One of the goals of the space-based gravitational-wave detector LISA [|l| is to make measurements that will probe 
gravity in the very strong field. The most stringent probes will involve binary black hole systems. Mergers of 
comparable mass black holes will be detectable throughout most of the universe j^; their measurements will probe 
the violent dynamics of two black holes combining to form a single hole. LISA will measure hundreds to thousands 
of wave cycles in such events, probing their waveforms with moderate precision. When the mass ratio of the system 
is extreme (M1/M2 ~ 10^'^ — 10^^), measurement of the waves probes the quiescent structure of the Kerr black hole 
spacetime. LISA should measure tens of thousands to millions of cycles from such inspirals (depending mostly on 
mass ratio) , and thus may be able to measure their waves — and probe black hole spacetimes — with extremely high 
precision. Extreme mass ratio inspirals will be the focus of this analysis. 

The case for studying extreme mass ratio inspirals has been argued at length elsewhere; I refer the reader to the 
Introductions of Refs. |^ and |^ for detailed discussion. Briefly, detailed studies have shown that compact stellar 
remnants in the central cusp of galaxies are scattered into very tight, eccentric orbits of the galaxies' black hole at a 
rate of at least several events per year out to a Gigaparsec (perhaps as high as an event per month). The final year 
or so of inspiral (depending on the large black hole's mass and spin, and on the binary's mass ratio) will occur as 
the small body spirals through the most extreme strong-field region of the black hole spacetime, from several horizon 
radii down to orbits just outside the hole's event horizon. 

By tracking the phase evolution of the waves over the many cycles radiated in this year, it should be possible to 
measure strong-field characteristics of the black hole spacetime. Fintan Ryan |^ has shown that these waves can be 
used to reconstruct the multipolar structure of the black hole. By the "no hair" theorem, a Kerr black hole's multipole 
moments are uniquely determined by its mass ill and spin a. Measuring more than two moments tests whether the 
massive object at the core of a galaxy is a black hole, or whether it is some more exotic compact object, such as a 
soliton star. 

Analyzing how well one can measure the multipoles of a massive compact object in this manner means computing 
gravitational waves from the inspiral of realistic small bodies on inclined, eccentric trajectories about bodies with 
arbitrary multipole moments. We are rather far from being able to analyze such a system. For now, I require that 
the massive body be a black hole — this greatly simplifies the description of its orbits. I also take the inspiraling 
body to be pointlike and nonspinning, sidestepping issues of spin-spin and spin-orbit coupling which complicate the 
inspiral and may in fact lead to chaos |9[[lCt] . Finally, because rigorous strong-field radiation reaction is currently being 
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developed , I use a "poor man's radiation reaction" approach, balancing the flux of energy and angular momentum 
in gravitational waves with the change in orbital energy and angular momentum. This formalism works well if the 
inspiral is adiabatic: the change in orb ital p arameters due to wave emission must be slow. (A more precise definition 
of adiabatic and "slow" is given in Sec. II B .) At present, the flux-balance approach cannot be used to flnd the change 
in the third constant of Kerr orbits (the Carter constant Q), so I require that the orbit be initially "circular" (of 
constant Boyer-Lindquist radius). These orbits do not become eccentric under adiabatic radiation reaction [|l2[-jl4[, 
constraining the system enough that Q's evolution can be inferred from the changes in the orbital energy and angular 
momentum alone. 

Given the large number of restrictions imposed, it is clear this analysis is rather far from the ultimate goal of laying 
the foundations for compact body multipole measurement with LISA. However, it is sufficient to begin developing 
gravitational waveforms and probing the inspiral trajectories for non-trivial orbits. These orbits are colored rather 
strongly by the orbital frequencies and their harmonics, and, despite the numerous restrictions, have a fairly ornate 
structure. This analysis is a useful way to produce gravitational waveforms that have a character similar to what we 
are likely to observe with LISA. They are restricted enough that they can be calculated with tools available now, but 
have enough structure that they can be used for developing data analysis tools and understanding how measurement 
of extreme mass ratio inspirals with LISA will work in practice. 

In the remainder of the Introduction, I summarize the structure and results of the paper. I first review the properties 
of circular, inclined orbits of Kerr black holes and the radiation reaction formalism used here. This material is presented 
in far greater detail in Refs. Q and |l6| . This paper just presents equations and points the reader to the relevant 
literature for a detailed derivation. In the end, radiation reaction data are presented as vectors (r, i) living at orbital 
coordinates (r, t) . These vectors describe how radiation drives the body to spiral through a sequence of orbits into 
the hole. The wave itself is developed from a set of complex amplitudes Zj^j^ which constitute a multipole expansion 
of the gravitational-wave strain measured by distant observers. This fully sets up this analysis. 

Inspiral trajectories are constructed by beginning at some coordinate (ro,io) and stepping inward along the data 
(f, 1) until the body hits the last dynamically stable orbit and plunges into the hole. The data are actually computed 
on a grid that covers the paramete r spac e of allowed orbits deep in the strong field, and are interpolated to get 
radiative effec ts off the grid. Section [III A| first describes how the data grid is set up, as well as which data are stored 
on it. Section [II B then briefly describes how the integration is done. 

Results are given in Sees. IV and First I discuss the inspiral trajectories in Sec. IV; they and their various 
features are displayed for several choices of black hole spin in Figs. |[ ||, P|, ^, and 0. During inspiral lasting 1-2 years, 
the inspiraling body executes several hundred thousand orbits in the extreme strong field of the massive black hole. 
This is not a surprise, and has of course been known for quite some time. More interestingly, the inclination angle 
barely changes during inspiral, particularly when spin a < 0.8M. This gives further credence to a suggestion by Curt 
Cutler that it would be worthwhile to explore "faking" the inspiral of generic Kerr orbits by holding the inclination 
angle fixed — this would fix the change in the Carter constant provided the changes in E and Lz were known. 

The influence of the black hole's event horizon on the inspiral is very strong. A very useful way to understand 
this effect is that the inspiraling body tidally distorts the black hole, causing the horizon to bulge. This bulge then 
exerts a torque back on the body. This viewpoint is explained at a very pedagogically accessible level in Chapter VII 
of Ref. ||r^ , which in turn is based on the work of Hawking and Hartle and Hartle jl^ . If the hole is rapidly 
rotating, this torque tends to significantly slow the inspiral, transferring some of the hole's rotational kinetic energy 
to the orbital motion. Inspiral can be prolonged for several weeks (~ 5% of the total inspiral time) by this effect, 
adding tens of thousands of additional orbits. By contrast, if the hole is rotating slowly, the torque exerted by the 
horizon bulge speeds up the inspiral. This fascinating coupling of the horizon to the inspiraling body's dynamics is 
an example of how strong-field features of the Kerr spacetime can be seen in gravitational waves. 

Section ^ discusses the waves produced as the body spirals inward. It should be strongly emphasized that computing 
waveforms is more difficult than computing inspiral trajectories. The waveforms require both the magnitude and phase 
of the complex amplitudes Zf^^j^ to high precision; the trajectories require just the magnitude. As such, the results 
shown in Sec. ^ should be taken as an important first step in generating extreme mass ratio inspiral waves, but can 
— and should! — be improved dramatically. 

Multiple harmonics of the orbital frequencies and fig are very important; cf. Fig. |^. Following a suggestion of 
Sam Finn, I write the waveform as a sum over many "voices" : h{t) ~ ^ hi„ik{t). Each voice hi„ikit) corresponds to a 
particular harmonic of the orbital frequencies, and has its own amplitude and phase evolution. Writing the waveform 
in this manner emphasizes that the phase evolution of each voice is rather simple, though the sum may be complicated 
and difficult to follow. It is likely this multi-voice structure will carry over to the generic case, adding an additional 
index for harmonics of the radial frequency fir- The analysis of LISA data may be facilitated by searching for extreme 
mass ratio inspirals voice-by-voice, rather than searching for the "chorus" of all voices at the same time. Concluding 



discussion, including suggestions for future work and extensions to this analysis, is given in Sec. Vl 



Throughout this paper, an overdot denotes d/dt, and a prime denotes d/dr. An overbar indicates complex conju- 
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gation. The quantities t, r, 9, and 4> refer to the Boyer-Lindquist coordinates. I use units with G = \ = c throughout. 



II. REVIEW: CIRCULAR ORBITS AND RADIATION REACTION 

A. Geodesic orbits 

Orbits of Kerr black holes are specified by choosing their energy i?, axial angular momentum L^, and Carter 
constant Q. The values used here have been divided by the small body's mass ^ [E and L^) or (Q) and are thus 
the specific energy, angular momentum and Carter constant. Having chosen these constants (and initial conditions), 
the orbit is then governed by geodesic equations for the small body's Boyer-Lindquist coordinates (t, r, 9, (f>) as a 
function of proper time r measured along its worldline; see ||20| . 

To describe circular orbits, it is useful to introduce the function R = Y?{dr/dT)'^: 

R = [E{r^ + a^) - aL,] ^ - A [r^ + (L, - aEf + Q] . (2.1) 

(The function Tj — -\- a? cos^ 6*, and A — 2Mr + a? .) Circular orbits satisfy i? = and R' — 0; stable circular 
orbits satisfy in addition R" < 0. Given the radius r and one of the three constants E, L^, or Q, it is straightforward 
to solve the system R = 0, R' = for the other two. In this paper, I choose r and E and then solve for and Q: 



, EM(r^ - a^) - AJr^(E^ - l)+rM , , 

\(a^ + r^)E - aLAr,E)] „ „ , , . 

Qir, E) = ^ ^^^^ - [r2 + a'E^ - 2aEL,{r, E) + L,(r, Ef] . (2.3) 

[Note that in Rcf. [|l6|, the a^E^ inside the second set of square brackets is incorrectly written a^E.] 

At a given radius, two orbits bound the behavior of all stable circular orbits. The most-bound orbit is the prograde 
equatorial orbit (t = 0). Its constants are [|2l[ 

E-^^ = ^-^-' + '1-' , (2.4) 
v/1 - 3^2 + 2qv^ 

^yl - 3i;2 + 2qv^ 

Q"* = , (2.6) 

where v = \/ Mjr and q = a/M . Prograde equatorial orbits exist at all radii outside Tpio, where 

rp„/M = 3 + Z2 + [(3 - Zi)(3 + Zi + 2Z2f'^ , (2.7) 

Zi = 1 + [1 - (a/Af)2] [(1 + a/M)i/3 + (1 - a/M)!/^] , (2.8) 

Z2 ^ [Sia/M^ + Zf]'^\ (2.9) 
This is also the radius of the LSO when t = 0° — no stable circular orbits exist inside rpio- If r < riot, where 

rret/A/ = 3 + ^2 + [(3 - Zi){3 + Zi + 2Z2)]'/' , (2.10) 

then the least-bound orbit is on the LSO, and the constants {E^^, L^^, Q^^) are found by numerically solving the system 

of equations R = Q, R' = Q, R" = 0. If r > Tret, the least-bound orbit is the retrograde equatorial orbit (t — 180°), 
and its constants are 

IK 1 — 2v'^ — qv^ 

S = = , 2.11 
v/1 - 3w2 - 2(7w3 

= -r.l±2£±£^ , (2.12) 
Vl - iv^ - 2qv^ 

Q"^ = . (2.13) 
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Orbits with Q 7^ are inclined with respect to the equatorial plane; a useful definition of the inclination angle l isQ 

cos 6 = . (2.14) 

Mapping out all stable circular or bits at some radius is fairly straightforward. First, calculate the constants that 



describe the most-bound orbit [Eqs. (|2.4D-(p.6|)]. Seco nd, c al culat e the constants for the least-bound orbit [solving the 
system i? = 0, i?' = 0, R" = if r < r^ot; using Eqs. ( p.ll )~(2.13) otherwise]. Finally, vary E from its most-bound to 



its least-bound extreme. 

Circular orbits are periodic, with two (generally incommensurate) frequencies fig and fi^, related to the small body's 
motion in the 9 and coordinates. These two frequencies (and their harmonics) strongly stamp the gravitational 
waveform. Typically, > fig; the differences are due to the oblate geometry of a rotating black hole and frame 
dragging (which can greatly augment fi^). A detailed derivation of these frequencies is given in Sec. lie of Ref. §. 

B. Radiative corrections 

The radiation reaction formalism used here is based on the Teukolsky equation 0, which governs the evolution 
of the complex Weyl curvature scalar related to radiative perturbations, '04. This section briefly outlines how the 
radiative corrections are computed, particularly the quantities relevant to this analysis. Reference Q discusses this 
material in greater detail. 

The "master equation" for the evolution of V'4 is separated with the multipolar decomposition Q 

1 r°° 

= ir-^acosBY E -.^.^ W^^^"^"-^ • (2-15) 

The t and dependences are trivial. The function -2Si^{6) is a spin- weighted spheroidal harmonic. It is useful for 
describing the 9 dependence of a radiation field with spin weight —2 in an oblate geometry. An effective algorithm 
for calculating this function is given in Appendix A of Ref. |^ . 

Computing the radial function Rimuj takes some effort. This function obeys the Teukolsky equation: 

2 d f 1 dRlmuj 



This equation is in self-adjoint form, and so can be solved with Green's functions [p3| . The homogeneous Teukolsky 
equation has two independent solutions, Rffni^ir) (which obeys the boundary condition that radiation must be purely 
ingoing at the event horizon) and R^^^ir) (which obeys the condition that radiation must be purely outgoing at 
infinity). From these solutions and from the source function Timuj{r), one can write down the general solution 

Rirm^ir): 

RimUr) = Z^,,^RZ.^{r) + ZZ^Rg,^{r) . (2.17) 

A detailed description of how one computes Zf^^^, Z^^, Rf^^, and R^^ (and how they relate to the source function 
%.mu}) is given in Ref. The most important details for this paper are as follows. First, for circular orbits, the 
source and all quantities derived from it arc describable as harmonics of fig and fJ^. Defining Lo^k — mfi^ + kVLg, the 
functions Zfj^^ and Z^^^ can be decomposed into harmonics: 

Z"r;::=Y.Z"r;::K^~^n.k). (2.18) 
k 



Second, as r ?-+ = M + %/ AP ~ (the coordinate of the event horizon), Zf^^^ and Z^^ constant; as 
r — > 00, Z^^ — > and Z^^^ constant. The behavior of the radiation at infinity thus depends only on Z^^j^, and 
at the horizon only on Z^f, . The gravitational waveform is 



^As discussed in Ref. this angle does not necessarily accord with intuitive notions of inclination angle. For example, except 
when a = 0, t is not the angle at which most observers would see the small body cross the equatorial plane. 
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thy, = J2 -^_25rm"' We''^"*'"'^'"'"*^ , (2.19) 



Imk ^'"/^ 



and the fluxes are 



dE 
dt 



ooM 




1 7^,00 12 
\^lmk 1 


rad 


\ ^ oo.H 

~ 2^ ^Imk 
Irak 


ooM 




1 ryH,00\2 


rad 


\ ^ oo,H 
~ ^Imk 
Irak 



(2.20) 



By assumption, the change in the orbital energy and angular momentum is opposite in sign to that radiated: E°'^^^^ = 
-E'^^'^, L°'''"* = —U^'^. Note the somewhat confusing reversal of H and oo on the Z/„ifc coefficients; this follows from 



the definition in Eq. (2.17). The factor ^ is a rather messy coefficient that follows from transforming the Kinnersley 
null tetrad (which is used to construct -ipi p^]) to the Hawking-Hartle null tetrad (which is well behaved on the 
event horizon); see Ref. for further details. (The factor Q!°° = 1; it is introduced so that these formulas can be 
written compactly.) As noted in there is a symmetry relation between quantities at {I, m, k) and {I, — m, —k): 

z"^::^k = • (2.21) 

This relationship can be used to reduce computation time — compute half the multipoles, use symmetry to get 
the other half. Alternatively, it can be used to improve computational accuracy — compute all the multipoles, use 
symmetry to reduce error: 



Imk 2 



^H,oo r^y+k^H.oo ^2.22) 

trak: comD ^ i — m~k: coniD V / 



This second approach was used here. 

A large number of terms were included when evaluating Eq. (2.20). A detailed discussion of the truncation criterion 
used here is given in Sec. VA of Ref. Q; in the language of that paper, — 10^** and e; = 10^''. The truncation error 
in E^'°° and L^'°° is therefore roughly 10"'^. Obtaining this level of accuracy required summing to at least I — 8 
(for innermost orbits of the a = 0.998A/ black hole, the sums were taken to I — 19). For each value of I, between 10 
and 30 values of k were included. In all cases, the sum over m ranged from —I to I. 

Using the Teukolsky equation to compute E^'°° and L^'°° requires that the inspiral be adiabatic: over orbital 
timescales, the trajectory must be nearly geodesic. This requirement enters through the source term, which depends 
on the inspiraling body's worldline. This is potentially the setup for a Catch-22: we need the worldline in order 
to compute the body's radiative corrections, but the corrections determine that worldline! The solution is to take 
the worldline to be geodesic for computing the source. We thus use the zeroth order (in n/M) geodesic motion to 
compute the first order radiative corrections. This is a good approximation only when the radiative change in any 
orbital quantity x over an orbit is much less than x- X^orb ^ X- Very useful quantities to monitor are the orbital 
frequencies, fl^ and fig. Because Tg = 2Tr/ftg, and defining = 2TT/n^, the adiabaticity condition for these quantities 
can be written 

1 ^Ig 

U^^,^ J^»l, (2.23) 
27^ n^y.g 

These parameters are closely related to the number of accumulated orbits in cj) and 6: 

N^,g J dti^4'.e = J dlnfl^^gAf^^g . (2.24) 

The condition A/^^e ^ 1 tells us that the orbits accumulated by the body as it passes through the orbital frequency 
band centered on fJ^^e and of width 5fl^^g = O^^e must be very large. More simply, the body must spend many orbits 
near any point in its orbital phase space. 

The Teukolsky equation can only tell us about the changes in E and L^- To fully describe the small body's spiral 
in, we impose circularity: adiabatic radiation reaction changes an orbit's radius and inclination angle, but it does 
not make the orbit eccentric JT^-Jl^ — circular orbits remain circular. Imposing this rule gives a relatively simple 
relationship between {Q,f) and {E,Lz); see Ref. |P, Sec. IIIA for details. Once Q and are known, it is simple to 
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get the rate of change of l using Eq. ( ^.14| ); cf. Eqs. (3.7) and (3.8) of Ref. ||^. The data (r, L) are the foundation for 
all the inspiral trajectories discussed in the next section. 

Letting x stand for any of the quantities [E, L^, Q, r), the total change in x is found by summing the change due 
to radiation flux out to infinity and due to radiation flux down the event horizon: 

x^iT+vx"- (2.25) 

The parameter rj (which clearly should equal 1) is introduced so that we can "turn off" the flux down the horizon. 
As we shall see later, setting 77 = is a very interesting test of how strongly the black hole's event horizon influences 
inspiral. 



III. COMPUTING INSPIRAL TRAJECTORIES AND WAVEFORMS 

Using the techniques described in Sec. ||, it is straightforward to compute the inspiral trajectory of a body orbiting 
a black hole and the gravitational waveforms generated in that inspiral. This is done in two steps. First, a "grid" of 
radiation reaction data is built in the phase space of allowed orbits in the black hole's strong field. The coordinates 
in this phase space are (r, /,). The radiation reaction data are the vectors along which gravitational- wave emission 
drives the orbit, (r,t), plus the a mplitude and frequency of the gravitational waveform, Zf^^. and LOmk- I describe 



how this grid is built in Sec. [II A. Next, the radiation reaction data are integrated to compute the inspiraling body's 



trajectory. The data are interpolated from discrete orbital coordinates on the grid to arbitrary points in the orbital 



phase space using a two dimensional cubic spline l22|. I briefly describe how this is done in Sec. HI B 



A. Making a radiation reaction grid 



Making a grid of radiation reaction data is for the most part straightforward. One discretizes the orbital phase 
space, choosing indices i and j such that radiation reaction data live at points (r,j , iij). Then, one uses the techniques 
summarized in Sec. || to compute the data. I describe here the algorithm used to discretize the orbital phase space, 
and also the actual form of the data (which are somewhat massaged so that they can be interpolated with as little 
error as possible). 

After some experimentation, I have found it convenient to make the grid evenly spaced in radius, and evenly spaced 
at each radius in the orbital energy. The discretized radius is thus 



J Sr. (3.1) 



The parameter rmin is chosen to be as close as is convenient to rpm [cf. Eq. (2.7)]. This creates a small gap in data 
coverage near the LSO at small values of l. I have used Sr — O.lAf in all computations; as will be discussed shortly, 
more intelligent choices can be made. The index j can in principle be made arbitrarily large. For the results discussed 
in Sees. IV and I set jmax ~ 20 — 30. For a 10^ M© black hole and an inspiraling body with /x = 1 Mq, this leads 
to an inspiral lasting about 630 days — typical observation time for a LISA mission. 
The discretized energy is 

E^j = £;,„b + l)SEj , (3.2) 

where SEj = [£'ib(rj) — E^\,{rj)] /(jmax ~ !)• I have used «max = 9. The parameter space coordinates {rj,Eij) fully 
determine the orbit at grid point (i,j). Using the relations given in Sec. Q, we can remap this coordinate to any 
convenient parameterization. I will typically write the grid points as {rj,Lij), with the understanding that this is a 
remapping from the coordinates {rj , Eij ) that are directly generated. 

A better choice for Sr can be made by studying how Z^^^i^ is computed from solutions of the Teukolsky equation. 
From Eq. (4.8) of Ref. g, we see that 

Zli^ (X r dr' Rg,^{r')%^Ur')/Mrr • (3-3) 



Here, Rf^^ is the solution to the source-free Teukolsky equation that is purely ingoing at the event horizon, Tim^{r) 
is the Teukolsky source term, and A(r) = — 2Mr -I- . The coefficient Zf^^, is then built from a harmonic 
decomposition of Zf^^. Both Timui and A vary relatively slowly with respect to r. The function Rf^^, on the other 
hand, oscillates with a phase factor that is roughly e"*^'' , where 
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2Mr+ r-r+ 2Afr_ r - r_ 
r (r) = r H in — — m — — — , (3.4) 

and where p = to — m^ln, the mode frequency modified by the hole's spin frequency = al2Mr+ [of. Eq. (4.4) of 
Ref. 1^]. The function r*{r) is the Kerr "tortoise coordinate", which often appears in studies of radiation propagation 
near black holes. Notice that dr* /dr gets quite large near the black hole. This tells us that e"'^'' begins oscillating 
very rapidly as the horizon is approached. The phase of Zf^^. is likely to change by ~ tt radians over a lengthscale 
5r* 1/pmk- This suggests laying out the grid evenly spaced in r* , or using a denser grid in r, with spacing 



5r 



dr* 



dr 



^ (3.5) 



\plk\ 



Here L is the maximum I index included in the calculation, and K is the maximum k index. This would create a grid 
that is far more densely sampled than that discussed here, especially for rapidly spinning holes. The rather crude 
choice Sr = O.IM places strict limits on waveform accuracy; this is an obvious starting point for improvements to this 
analysis. 

Some care must be taken to put data on the grid that are useful for generating the inspiral trajectory. For the 
radiation reaction quantities, I have found it useful to store d cos t/dt and 

[cosL - cosLi^so{r)]r = p , (3.6) 

rather than (f, L). The quantity d cos o/dt is useful simply because it is more naturally related to the rates of change 
and Q. The quantity p on the other hand accounts for the fact that f diverges as the LSO is approached — 
the prefactor cosi — costLSo('') nicely clears out this divergence. This makes it possible for r to be very accurately 
interpolated from i = 0° to tLSo(^)- For the gravitational waveform, the frequencies fi^ and Qg are stored, as are the 
magnitude Aimk and phase ^imk of the complex amplitude, 

Zf^^i. = Ainik exp(i$/„fe) . (3.7) 

For equatorial orbits, Zl^^, = when k 0. The numerical code which computes ^imk from ^j^j, erroneously assigns 
^imk = in this case. This is because the code cannot evaluate ^imk — arctan[Im(Z;^j.)/Re(Z;^jj,)] in the hmit 
^imk ~* ^- "Fo get around this problem, ^i^m,k^o is only stored for non-equatorial orbits. Near the equator, ^i^m,k^o 
can then be computed quite well using extrapolation. 

An example grid, for a hole with a = 0.998A/, is shown in Fig. |l]. The arrows in this plot represent the vector (f, L). 
Notice that the evolution gets more rapid as the LSO (represented by the dotted line) is approached — the arrows 
get longer, until finally the orbit becomes dynamically unstable and plunges into the hole. (The diverging vectors on 
the LSO have been suppressed in the figure, as they tend to overwhelm the rest of the data.) This grid prototypes all 
the radiation reaction data used in this paper — although differing in detail depending on the black hole spin, they 
have the general shape and layout shown in Fig. |^. 

B. Integrating a trajectory across the grid 

Radiation reaction data on the grid are used to build inspiral trajectories — the paths [r{t), t(t)] that bodies follow 
as they spiral into the black hole. This is done with Euler's method: given a starting point {tq, lq] and a time step 5t, 
advance to [ro + rSt, lo + iSt]. It would be straightforward to use a more sophisticated stepping algorithm, such as 
a Runge-Kutta method, to improve the accuracy of the inspiral trajectory. For the purpose of a first exploration of 
extreme mass ratio inspiral properties, Euler's method is a good compromise between accuracy on the one hand, and 
code complexity and computational cost on the other. As the trajectory is built, other useful data are also calculated, 
particularly the accumulated number of orbits and the gravitational waveform. 

The Euler method integration requires derivatives at arbitrary points in the orbital parameter space, which in turn 
requires smooth methods for interpolating off of the grid points. Linear interpolations turn out not to be adequate: 
discontinuities in the data's second derivative as grid boundaries are crossed noticeably affects the phasing of the 
computed gravitational waveforrrQ Two dimensional cubic spline interpolation | |2^ seems to produce an inspiral 
trajectory that is adequately smooth; it is used for all interpolations in this analysis. 



^The discontinuities are amazingly clear when one transforms a gravitational waveform so produced into audio data. 
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The spline interpolation is done in the index coordinates, Suppose we wis h to interpolate data for some 

field X onto (r, t). First, we interpolate in j [which maps simply onto r; cf. Eq. ( |3.l| )] for each value of i, yielding a 
one-index set of data at the radius r: 

Spline in r : Xij ^ Xii^) ■ (3-8) 

Since all of the relevant data now live at a single radius, the index i maps onto cos t. Interpolate in cos t to get the 
final data: 

Spline in cost : Xiif) xi''^^'-) ■ (3-9) 



This procedure behaves very well with the radiation reaction data and waveform data discussed in Sec. [II A 



The procedure for generating the inspiral trajectory and gravitational waveform now reduces to a simple recipe: 

1. Pick a starting coordinate, (tq, LQ;t = 0). 

2. Interpolate the data dcosi/dt, [cost — cosiLSo] ^ = P^ ^(t>^ -Aimk, and $;„ifc onto this coordinate. Note that 
for the waveform amplitude and phase, this requires interpolations for every value of /, m, and fc, which becomes 



computationally intensive. As discussed in Sec. [II A, $;„ifc is not stored on the equator when k ^ 0. One gets 
data for ^L,n,k^o between the equator and the next grid point by extrapolation. This works very well because, 
at constant radius, $;mfc is nearly flat as a function of cos l. 

3. Generate the gravitational waveform at that moment on the trajectory: 

i^mk = + kfle , 

h{t) = h+it) - ihy. (t) J2 ^^-^^Z"'" ^^y'f^r..e^■m^e-2^^irnN^+kNe) (3^0) 

1=2 m=-l k=-k„-,^^ ^'n*: 

Here, d is the angle between the observer's line of sight and the hole's spin axis, and (p is the orbit's (f> coordinate 
at t = 0. Formally, Zmax = fcmax = oo. In practice, these numbers are truncated at something much smaller. The 
last term in h{t) is a harmonic of the accumulated orbital phases; the factor in the exponential is equivalent to 
—i J dt{mfl(j, + kQg). It would reduce to ~2Tri{mQ^ + kQ0)t if the orbital frequencies did not themselves evolve. 

4. Compute the number of orbits accumulated so far: 

N^t + St) = N^{t) + 5t/2iT , 

Ne{t + 5t) = Ng{t) + VLg5t/2T: . (3.11) 

Since the frequencies of and 9 motion are generally incommensurate, these two numbers can be quite different. 

5. Take an Euler step to a new coordinate: 

^ncw — ^old St , 

'"new = ?'old H P5t , 

cos i — COS iLSO 

dcOSL ^, /o -,f,^ 

COS tnow = COS told H ; — St. (3.12) 

at 

6. Go to Step 2 and repeat. Gontinue until the inspiraling body crosses the LSO. 

Results from applying this analysis to a large number of inspirals of black holes of various spins are discussed next. 



IV. RESULTS: INSPIRAL TRAJECTORIES 



Before generating inspiral trajectories and gravitational waveforms, we must choose a sample of black hole spins 
— the trajectories are unique for each spin choice of the massive black hole. (The effect of varying other parameters, 
particularly the mass and mass ratio, can be obtained by scaling.) We would like a range of a that covers at least 
qualitatively the range likely in astrophysical extreme mass ratio inspirals. 
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This analysis was done for four representative black hole spins: a = 0.998M, a = 0.3594M, a = 0.8M, and 
a = 0.05M. The first two were selected because they are rigorously calculable. The value a = 0.998M is the 
"astrophysically maximal" spin one finds when the black hole evolves due to thin disk accretion: preferential capture 
of counter-rotating photons versus co-rotating photons buffers the spin at 0.998M, preventing the hole from reaching 
the Kerr maximal value a — M [ p5| . The second choice is obtained by locking the rotation frequency of the horizon, 

= a/2A/r_|_, to the orbital frequency of the innermost equatorial orbit at r-pro, ^4, — M /[rpro + aAI^^^]. Equating 



these two frequencies and using Eq. (2/7) yields a — 0.3594M. This spin might be obtained if at some point in 
the black hole's history strong magnetic fields threaded the event horizon and the inner edges of an accretion disk, 
torquing the hole such that it became locked to the disk's rotation |2^. The third choice, a = 0.8M, is primarily 
chosen because it breaks up the large gap between 0.998A/ and 0.3594M. It is worth noting that this spin is in the 
range predicted by detailed models of black hole evolution in the presence of magnetohydrodynamic torque, such as 
are described in some models of quasar engines [2^. Finally, a = O.OSAf is chosen to give an example of inspiral 
into a slowly spinning black hole. In all cases, the numbers presented are for a /i = 1 Mq body spiraling into an 
M = 10^ Mq black hole. 

Results are summarized in Figs. H, |[ |5[ ^, and ^ Before discussing each case, it is worth noting some general trends 
in the data. The span of data in all cases has been chosen so that the total inspiral lasts about 630 — 650 days for 
the shallowest inclination angles. Inspiral typically halts when the small body crosses the LSO; at the very shallowest 
inclination angles, it halts when the span of radiation reaction data ends (very close to the LSO). 

The radius of the LSO as a function of inclination angle, rLso('^)i depends quite strongly on the black hole spin 
— consider that for a ~ 0, rLso('') = 6A/ independent of l, whereas for a — M, r-Lso(0°) — M, rLso(180°) — 9M. 
Because of this, the starting radius of the inspiral is moved out to larger values as a is decreased. Also, the total 
number of starting inclination angles is less for large a: when a is close to M , trajectories that begin at high values 
of L hit the LSO more quickly. Thus, in the data discussed below, inspiral into small a black holes occurs at larger 
coordinate radius than for large a, and more starting inclination angles are included in those data sets. 

In all cases, I show the accumulated inspiral time, the number of orbits about the spin axis {N^) and the accumulated 
number of oscillations in 6 (Ng). These numbers are easily scaled to different masses and mass ratios: the trajectory 
shapes are independent of the masses (provided the radial axis is r/M), and 

M 

Tinspiral OC , N^_g OC . (4.1) 

The trajectories shown here thus apply for any black hole mass and any extreme mass ratio. 



A. Trajectories for a = 0.998M 

Trajectories for inspiral into a hole with spin a = G.998M are summarized in Fig. |^. This figure shows inspiral 
from r = AM into the LSO for several starting values of the inclination angle l. The accumulated inspiral time, orbits 
about the spin axis, and number of oscillations in 9 label each trajectory. 

Notice that the trajectories are nearly flat — l barely changes, decreasing somewhat as the small body passes 
through the very strong field. This decrease is contrary to weak-field expectations djlj], and arises because the 
characteristics of Kerr geodesies become locked to the event horizon in the very strong field; see Ref. [|l^ for detailed 
discussion. The fact that the change in t is so small is very interesting. The change turns out to be even smaller when 
a is not so large. 

Figure || illustrates the true strong-field inspiral sequences predicted by general relativity, including the effects of 
radiation out to infinity and down the event horiz on. An interesting experiment is to "turn off" the radiation down 
the horizon, setting the parameter 77 in Eq. ( ^.25| ) to 0. The inspiral trajectories found in this exercise are plotted 



in Fig. ^. Aside from turning off the horizon flux, these trajectories are computed with conditions identical to the 
trajectories plotted in Fig. ||. Notice that, particularly for shallow inclination angle (6 ^ 50°), the small body inspirals 
more quickly when the flux down the horizon is not included. The effect is quite significant: inspiral of the shallowest 
orbits is shortened by several weeks, orbiting ^ 10^ fewer times when the horizon fiux is ignored. 

That the horizon fiux tends to prolong the inspiral is, on first consideration, extremely surprising. Intuitively, one 
would guess that each flux would be an energy sink, so that each reduces the orbit's energy. Because of the hole's 
rapid rotation, this is not the case: the phenomenon of superradiance plays a major role. Superradiance is essentially 
a manifestation of the Penrose process — a radiation pulse incident on the black hole is backscattered with a gain 
in the pulse's energy |Q. The energy gain is at the expense of the hole's rotation — rotational kinetic energy is 
converted to radiative energy. Thus, when the hole spins rapidly enough, the horizon is a source of energy, not a sink. 
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A rather more physical picture of this phenomenon (which explains why the energy is transferred to the orbit, not 
simply radiated to infinity) can be put together based on work by James Hartle . The portion of the Weyl tensor 
which describes radiation down the event horizon can be understood as the tidal field of the orbiting body acting on 
the black hole. This field distorts the horizon, raising a tidal bulge on the hole. The distortion can be quantified by 
computing the curvature of the horizon. Using formulas from Rcf. jl^ , it is easy to show that the distortion is given 
by a sum over the coefficients Zj^^^. that in this analysis describe the radiation flux at the horizon. 

The black hole's tidal bulge is equivalent to the bulge that the moon raises on the earth. When the rotation 
frequency of the hole is greater than the frequency, Qh > ^4>, the hole's spin drags the bulge ahead of the orbiting 
body, in exactly the manner in which a tidal bulge on a viscous, fluid body is dragged by the body's rotation. (One 
can in fact interpret the horizon's behavior in terms of a viscous membrane |l7|| .) In a reference frame that co-rotates 
with the small body's orbit, the bulge appears to lead the small body by some angle. The bulge thus exerts a torque 
that acts to increase the body's orbital velocity, transferring the hole's rotational energy into the orbit and slowing 
the inspiral. A very accessible discussion of this point (which directly relates the "tidal bulge" effect to a torque) is 
given in Chapter VII o f Ref . pTf ; the fundamental underpinnings of this analysis are developed in Refs. |T^jl9| ]. 

As discussed in Sec. II B , it is important to monitor the adiabaticity parameters A/",/, and Afe during the inspiral. 
The values of A/^ and A/g for the inspiral beginning at i = 20° are plotted in Fig. |4| A/^ is plotted as the solid line, A/g 
is the dotted line. We see the adiabaticity condition A/^.e ^ 1 is strongly satisfied everywhere except right before the 
end of inspiral. (This is because the inspiral rate sweeps up as inspiral proceeds, and can be regarded as a precursor 
to the final plunge.) The adiabaticity parameters scale with M/^. This is a warning that the results of this code 
cannot be taken seriously if the mass ratio is not extreme enough: when M/fi < 1000, the final stages of inspiral 
will not be adiabatic. An interesting feature in this plot is the divergence in Afg near t = 470 days. This divergence 
occurs because fig switches sign there: the frequency stops increasing and begins decreasing. This will be discussed 
at greater length in Sec. |^. 

The adiabaticity parameters appear qualitatively as in Fig. I in all cases studied here (with the modification that 
the divergence in Afg only occurs for a = 0.998M), so no other examples of A/^,© will be shown. 



B. Trajectories for a = 0.3594M 

Trajectories for inspiral into a black hole with spin a — 0.3594^ are summarized in Fig. ^. These trajectories are 
set up in a similar manner as those for a = 0.998M; notable differences are that they begin dX r — 6.6M, and more 
starting values of i are included. The trajectories are even more fiat in this case than they are when a = 0.998M. The 
inclination angle actually increases slightly in all cases (as weak-field analyses predict), but the amount of increase is 
too small to be noticeable on the figure. 

Turning off the horizon flux does not have a very marked effect on inspiral for this spin, as can be seen by 
comparing the upper panel (horizon flux included) and lower panel (horizon flux ignored) of Fig. ^ This is essentially 
by construction: the spin a = 0.3594M is the value for which the horizon spin frequency matches the innermost orbital 
frequency. The lead angle between the tidal bulge and the orbiting body ranges from very small to nonexistent — as 
the body spirals through the innermost orbits, the bulge is nearly perfectly lined up with the body, exerting almost 
no torque on it. The small difference between the trajectories in the two panels is not surprising. 



C. Trajectories for a = 0.8M 

Inspiral trajectories for a black hole with a = O.SAf are shown in Fig. ||. Not surprisingly, the properties of these 
trajectories are intermediate to those shown for the cases a = 0.998M and a — 0.3594M. As in the case a — 0.3594M, 
the change in inclination angle is extremely small, increasing slightly. 

Inspiral is faster when the horizon flux is not included. This, again, is not a surprise: we expect that when a = 0.8M 
the bulge raised on the event horizon leads the orbiting body, so that tidal coupling transfers energy from the hole's 
spin to the orbit. Although qualitatively the effect is the same as for the case a — 0.998Af , quantitatively the effect 
is much smaller. The fractional change in the inspiral time is no more than about 1% when a = 0.8M, as opposed 
to 4.2% when a = 0.998A/. Part of this is simply because the spin frequency of the hole is smaller when a — 0.8M 
— the tidal bulge does not lead the orbiting body quite as much. Also the body's orbits never come as close to the 
event horizon as they do when a = 0.998M — they reach the dynamical instability before they get so close. As a 
consequence, the tidal coupling is never as strong, and the integrated effect on the trajectory is relatively small. 
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D. Trajectories for a = 0.05M 



Figure |^ shows trajectories for inspiral into a hole with a = 0.05A/. The dependence of the inspiral properties 
on inchnation angle is weak compared to the other cases examined, though definitely present. The change in the 
inclination angle is practically unnoticeable; the largest change (for inspiral near 90°) is Si ~ 0.007°. 

In this case, inspiral is quicker when the horizon flux is included. The horizon functions as a sink of energy, in 
accord with simple intuition: when the horizon flux is not included, the small body does not spiral in so quickly. 
Because the hole rotates so slowly, the tidal bulge raised on the horizon lags the orbiting body, and the torque that 
it exerts on the orbit tends to increase its inspiral rate. The magnitude of the effect is quite small (fractional change 
in inspiral time is about 0.1%). 



V. RESULTS: GRAVITATIONAL WAVEFORMS 



As the inspiral trajectory is generated, the gravitational waveform is constructed using Eq. ( 3.10| ). An example 



waveform is shown in Fig. |8|. This is the + polarization generated during the inspiral beginning at t — 40° into a 
hole with a — 0.998Af; cf. Fig. ||. It is viewed in the hole's equatorial plane {9 = 7r/2). The mass ratio used here is 
10~^; recall that 10~^ was used to make Fig. ^. The increase was to speed up the inspiral so that the amount of data 
generated in the wave was kept manageable. The waveform contains contributions from ^ = 2, 3, 4, and k = —4 to 4. 

Some of the interesting features of the waveform are apparent in Fig. ||. The modulation of the "carrier" signal by 
orbital motion in is plain, as is the evolution of the signal frequencies. Some of these features are made even clearer 
by transforming the gravitational waveform into audio data and playing the gravitational-wave "sound" . I have placed 
sounds corresponding to such gravitational- wave signals at the URL listed as Ref. Q ; the reader is invited to play 
these sounds and judge for themselves how clear are the various inspiral features. (In order for the sounds to lie 
within the frequency band to which human ears are sensitive — and to keep the signal duration reasonable — the 
mass of the black hole is scaled down to M ~ 100 Af©. Discussion of this and other technical points is given on the 
page listed in Q.) 

Extreme mass ratio gravitational waveforms are very usefully described as multi-voice chirps. Strong radiation is 
emitted at multiple harmonics of the inspiral frequencies; each harmonic plays the role of a separate "voice" in the 
"chorus" that is the overall waveform. This structure is quite apparent in the audio versions of the waveforms. It can 
be made clearer by rewriting the waveform Eq. ( ^.ig| ) as follows: 

'max I fcniiix 

h{t;^,^) = J2 J2 7^^™fe(^;^?)exp[^^'^™fc(^;^)] . (5.1) 

1 — 2 rn— — l k— — k^-^x 



Comparing with Eq. ( S.lOj ), we read off the values of the wave amplitude Tiimk and phase ^'jmfe in terms of quantities 



directly computable from the Teukolsky equation: 

Uir,At;^) - Ai^k-2SZ'-'{^)l^lk • (5.2) 

"^Imkit^f) ^ ^Imk + ITlLp - dt' {mVl^ + k^g) . (5.3) 







[Recall Eq. (3.7), defining Aimk and ^imk] The last line emphasizes that a big portion of the phase accumulated is 
nothing more than the integrated orbital phases. 

Presenting the waveform in this way gives a very good sense of which harmonics contribute strongly to a 
gravitational-wave measurement. As will be discussed below, the phase evolution for each voice is very simple, 
though their summed effect is typically rather complicated (cf . the summed waveform shown in Fig. |^) . This suggests 
that searching for each individual voice may be more effective than searching for the entire waveform. An effective 
way of implementing such a search might be developed using the "Fast Chirp Transform" (FCT) of Jenet and Prince 
p9[ . The FCT is a generalization of the Fast Fourier Transform (FFT; see, e.g., Ref. and references therein). The 
EFT provides an efhcient computational implementation of a Fourier transform — a decomposition of data onto basis 
functions e^^*-''*. The FCT generalizes this to decompose data onto functions whose phases do not vary linearly with 
respect to time. The phase behavior of the basis functions e**'''^'^^*'' must be specified by some model; indications are 
that it works well for a wide-variety of smoothly varying phase behaviors. Although further investigation is needed 
to test this idea, the voices of extreme mass ratio inspiral are good candidates for detection with the FCT. 

The results here are restricted to a = 0.998Af and a = 0.3594Af. The results for these spins contain enough 
interesting structure to draw some useful conclusions about extreme mass ratio gravitational waves. 
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A. Waveforms for a = 0.998M 



Discussion in this section wiU focus on waves generated during an inspiral that starts at t = 20° (cf. Fig. All 
data shown corresponds to waves measured in the hole's equatorial plane. 

Figure || shows the amplitudes Tiimk and phases "^imk for Z = 2, m = 2, and k g [0,1,2,3,4]. Notice that the 
accumulated phase 5'22fe for all values of k is quite large — the number of accumulated radians ranges from about 
8 X 10^ to 1.1 X 10^. The accumulated phase scales inversely with the mass ratio, so the accumulated phase in these 
harmonics would be ~ 10^ radians for a mass ratio of 10~^ (as was used in Figs. ||, ^, ^, |^). It is by tracking the phase 
over this large number of radians that it will be possible to determine the characteristics of the black hole spacetime 
with very high accuracy. 

The strongest waves are emitted in the fc = harmonic. Indeed, the I — 2, m — 2, k ^ voice turns out to be the 
strongest of the various harmonics that contribute to the gravitational waveform — not too surprising, since this is 
the leading quadrupole component of the waves. Note, though, that the k = 1 harmonic becomes almost as strong 
as the k = voice late in the inspiral. This is a generic feature of the multi-voice chirps: late in the inspiral, when 
waves are emitted from deep in the black hole's strong field, harmonics which were not initially very important can 
contribute strongly to the waves. 

The amplitudes corresponding to fc = 2, fc = 3 and fc = 4 "crash" and behave rather oddly late in the inspiral. 
This is a resolution problem caused by the crude manner with which I have discretized the orbital parameter space — 
constant Boyer-Lindquist radial st eps with 6r = Q.IM. In this region, pmk ~ (0.2 — 0.5)/M and \dr* /dr\ ~ 30 — 40. 
Following the discussion near Eq. (3^), the choice Sr = 0.1 A/ is not really sufficient to follow the phase variation of 
^imk'-i decreasing the step by a factor of 2 to 5 would vastly improve this analysis. 

The behavior for I = 2, m = is plotted in Fig. |l^. These phases actually decrease at late times. The strongest 
contributor to the phase evolution of these waves is the frequency ^Ig, which begins decreasing at late stages of the 
inspiral. The evolution of fig and fl^ are shown in the upper panel of Fig. |ll|. The decrease in ^Ig occurs because 
the small body orbits very close to the event horizon towards the end of inspiral. Its 6 motion is highly redshifted as 
seen by distant observers. (The (j) motion, by contrast, is not redshifted; instead, it limits to the hole's spin frequency 
as the orbit becomes locked to the dragging of inertial frames.) No reverse chirp occurs when a = 0.3594M — the 
inspiraling body is never close enough to the horizon for there to be a significant redshift effect. 

Measuring a reverse chirp would be in principle a wonderful probe of the black hole's strong field properties, since 
it can only occur for orbits that get very close to the event horizon. Unfortunately, the reversal in the phase evolution 
coincides with a rather rapid decrease in amplitude for these voices, as can be seen by comparing the upper and lower 



panels of Fig. 10. This is because the relevant source multipole moments vary more slowly as the frequency drops. 
Since this variation sets the wave amplitude, any voice with a decreasing frequency will likewise have a decreasing 
amplitude pOf . 

Although the waveform's strongest radiation is emitted in voices corresponding to I = 2, m = 2, the voices 
corresponding to Z = m, to ^ 2 are also quite strong, particularly at late times. For example, at the end of inspiral, 
the / = 4, TO = 4, fc = voice is about one fifth as strong as the I — 2, m — 2, k — voice. 



B. Waveforms for a = 0.3594M 

I again focus on waves generated during inspiral beginning at t = 20° (cf. Fig. The phase and amplitude 
evolution are qualitatively similar to those shown in Fig. ^ so I will not show additional plots. I will focus on voices 
for ; = 2, TO = 2, and k £ [0, 1,2,3,4]. 

As when a = 0.998Af , the voice corresponding to /c is strongest. Unlike that case, no other voices become very 
strong towards the end of inspiral — when the inspiral ends, the k = 1 voice is about a factor of 10 weaker than the 
k = voice. For this spin, the code was able to reliably compute all voices except the one corresponding to fc = 4. 
This is because the frequencies are quite a bit lower in this case, and \dr* /dr\ < 1.7 over the inspiral domain. The 
grid spacing Sr = 0.1 Af works well, except for large harmonic indices. 

Two important points should be noted from the results for other voices. Because fig grows monotonically, no voices 
fade away, in contrast to the Z = 2, to = voices for a = 0.998A/. Also, the voices for Z > 2 do not become as strong 
at late times as they do when a = 0.998Af . At late times, the Z = 4, to = 4, fc = voice is about one twentieth the 
strength of the I — 2, m = 2, k — voice. 
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VI. CONCLUSION: IMPLICATIONS FOR LISA SOURCES 



The results presented here give the first strong-field gravitational waveforms for rigorously computed extreme mass 
ratio inspiral trajectories. Because these trajectories are restricted to zero eccentricity, they do not correspond to 
inspirals that LISA is likely to observe. Nonetheless, they contain interesting features that are very likely to be present 
in the general case, and are certainly worthy of further study. 

One of the most interesting features of these extreme mass ratio inspiral waves is their "multi-voice" structure: 
each (/, m, k) harmonic follows its own phase and amplitude evolution. The total waveform is given by the sum of 
the various voices as the inspiral progresses. Provided inspiral is adiabatic, this multi-voice structure should describe 
generic inspiral waves as well, with a third index describing harmonics of the radial frequency fir- 

It should be possible to take advantage of the multi- voice structure of extreme mass ratio inspirals when developing 
strategies for analyzing the LISA datastream. Searches of data from ground-based detectors such as LIGO will rely in 
many cases on matched filtering, a technique that cross-correlates the instrumental data with templates of a source's 
gravitational waveform. Matched filtering works particularly well when the template accurately models the source's 
phase evolution; accurate models for the amplitude are far less important. In the case of extreme mass ratio inspiral, 
it may be far easier to search for each voice with its simple phase evolution, — ^^m/c(Os**''"'°*■*^ than to search 

for the entire "chorus", h(t) — himk{t), with its comparatively complicated phase evolution. Some signals might 
only be detectable in their higher harmonics — one can imagine using the m = 4 voices to find inspiral into holes 
so massive that the m = 2 harmonics are at frequencies too low to be seen. Techniques that look for extreme mass 
ratio inspirals on a voice-by-voice basis should have no problem finding those inspirals, provided these high-harmonic 
voices are not too weak. The Fast Chirp Transform [E9| may be a computationally eiTective means of implementing 
such a voice-by-voice search. 

In this vein, it is interesting to no te that the phase evolution of each voice is dominated by that voice's integrated 
frequency. Recall, from Eqs. ( ^.19| ) and (5T) ~ (x3) that the evolving phase consists of two pieces (ignoring the 



constant offset mip, which arises from initial conditions). One piece amounts to an integral of that voice's frequency, 
LUmk, and contributes ^ 10^ — 10^ radians to the inspiral. The other corresponds to the phase of the complex 
amplitude ^/^j, at a given moment; it only accumulates a few to a few tens of radians. Accurately calculating the 
inspiral trajectory, and hence having good information about the evolution of the orbital frequencies, is likely to be 
more important than knowing $/mfe. 

The influence of tidal coupling between the inspiraling body and the black hole is very robust and interesting — 
such coupling will strongly influence the inspiral. If nature provides black holes that rotate sufficiently fast, the 
prolonging of inspiral due to this coupling will easily be seen in the gravitational- wave data. This is a clean, beautiful 
probe of the strong-field black hole spacetime. 

The analysis presented here suggests several possible directions for future work: 

• Equatorial orbits. Glampedakis and Kennefick ||l5[| are currently finishing an analysis of radiation reaction on 
eccentric, equatorial orbits of Kerr black holes; it is an analysis equivalent to that done in Ref. ||^. The radiation 
reaction data they develop could be used to duplicate this analysis but to find inspiral sequences for eccentric, 
equatorial orbits and to develop the associated inspiral waveforms. It would be quite valuable to have results 
that include the effect of eccentricity and modulations induced by the radial frequency fir- 

Generic orbits. We can in principle calculate gravitational waves from generic orbits of Ker r black holes — all 



we need to do is specify the source term for a particular orbit in the Teukolsky equation (2.16) and solve. We are 
hampered, though by our inability to compute backreaction on such orbits. Valuable insight may be developed 
by imposing physically motivated constraints to produce an approximate inspiral trajectory. The trajectories 
presented here indicate that the change in inclination angle is very small. Using the Teukolsky equation to 
compute E and from the gravitational-wave flux and imposing i = is sufficient to fully compute the inspiral 
trajectory through the generic orbit parameter space (r, e, (.). The gravitational waves developed in such a study 
would be valuable tools for exploring the influence of all three orbital frequencies (Jl^, Og, fir) on the waveform, 
even if the waves so produced do not correspond exactly to those in nature. 



Improved waveforms. As discussed in Sec. Ill A, the parameter space discretization used here was rather crude. 
It was in fact insufficient to accurately compute the amplitudes and phases of high frequency contributions to 
the waveforms shown here. As work begins to develop data analysis tools for LISA, this analysis should be 
refined so that simple limitations such as this discretization are unimportant. Also, the parameter space covered 
in this analysis is rather inadequate for studying LISA measurements. The range covered here was selected 
simply for computational convenience. A better analysis would chose the parameter space covering to span a 
range of inspirals corresponding to the masses and observation times that LISA is likely to measure. 
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Implementing the above items will put on us well on the way to understanding the data analysis problem with 
LISA. Having a broader set of accurate waveforms will serve as a very useful testbed for developing analysis tools and 
for getting a better understanding of how measurable extreme mass ratio signals are likely to be. 
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[30] In several talks, I have claimed that the amplitude associated with the decreasing frequency will actually be rather large. 
This was entirely wrong. My statement was based on an earlier stage of this analysis; at the time, the code which interpolates 
from the parameter space grid to arbitrary parameter space coordinates had a significant bug, which introduced a strong, 
backward chirping component into the waveform. 
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FIG. 1. Circular orbit radiation reaction data near the last stable orbits (LSO) of a Kerr black hole with a = 0.998M. Any 
point (r, t) in this plot is a circular geodesic orbit. The arrows are proportional to the vector [{M/fi)r, (M'^ /fi)L]: the orientation 
gives the direction in which gravitational- wave emission drives the orbit, and the magnitude is proportional to the rate at which 
it is so driven. The dotted line is the LSO — orbits above and to the left of this line are dynamically unstable and rapidly 
plunge into the black hole. The arrows get longer as this line is approached and their stability decreases. (Additional data 
were produced representing radiation reaction for orbits on the LSO; these orbits are so unstable that their radiation reaction 
vectors do not fit on this plot. These data are used in all computations, however.) 
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FIG. 2. Inspiral trajectories in the strong field of a Kerr bl ack h ole with a = 0.998M. To make this plot, the data shown 
in Fig. |l|were integrated using the procedures discussed in Sec. [IIB , assuming that the black hole has mass M — 10^ Mq and 
that the inspiraling body has mass = 1 Mq . The trajectory shapes are independent of the two masses, so the inspiral times 
and accumulated number of cycles N^.g can be rescaled to other masses quite easily: Tinspirai oc M'^ / fj,, N^^g <x M/^. Notice 
that the trajectories are nearly flat — l decreases, but not very much, over these inspirals. 
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FIG. 3. Insp iral trajectories in the strong field of a Kerr black hole with a = 0.998M, ignoring the flux down the horizon 
— cf. Eq. ( ^.25[ ), with rj = 0. The figure is otherwise identical to Fig. ^ The trajectory shapes change very slightly (t does not 
decrease quite as much over the inspiral), though that effect is very small. Much more interestingly, the inspiral is markedly 
faster, especially at shallow inclination angle, the s mall body takes several weeks less spiraling to the LSO, and executes many 



thousands fewer orbits. As discussed in Sec. [V A , this illustrates how tidal coupling between the small body and the event 
horizon strongly impacts the inspiral: a tidal bulged is raised on the hole, which, due to the hole's rapid rotation in this case, 
transfers rotational kinetic energy to the small body's orbit. 
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FIG. 4. The evolution of the adiabaticity parameters Af,p^e ~ ^^^,s/27rf2^_e for the inspiral track beginiimg at t = 20° shown 
in Fig. 1^. The solid line is A/",^, the dotted line J\fe- Except at the extreme end, A/^.e ^ 1, indicating that the inspiral is indeed 
adiabatic, as required. Note that A/',^,e oc M/jj., indicating that the adiabatic requirements are unlikely to be met if M/fi < 1000. 
The divergence in Ae is due to Vlg passing through zero and changing sign during the inspiral — it ceases chirping up, and 
begins chirping down. 
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FIG. 5. Inspiral trajectories in the strong field of a Kerr black hole with a = 0.3594Af . The top panel includes the effects of 
the down-horizon flux; the bottom panel does not. The span of data is chosen so that the total inspiral duration is similar to 
that shown in Fig. ^. For this spin, inspiral is nearly identical with the horizon flux included or disincluded: inspiral is slightly 
faster without horizon flux, but not nearly so much faster as when a = 0.998Af. This is largely because the spin is not fast 
enough to drag the tidal bulge on the horizon so far forward. In fact, at the innermost orbits the bulge should be essentially in 
perfect face-on lock with the orbiting body, since a — 0.3594M is the spin value at which the horizon spin frequency matches 
the innermost orbital frequency. 
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FIG. 6. Inspiral trajectories in the strong field of a Kerr black hole with a — 0.8M. The top panel includes the effects of 
down- horizon flux; the bottom panel does not. The span of data is chosen so that the the total inspiral duration is similar to 
that shown in Fig. ^. As when a — 0.998M, inspiral is quicker when the horizon flux is ignored. However, the magnitude of 
the effect is much smaller. The torque on the orbit that the horizon's tidal bulge exerts is not as great because the hole rotates 
less quickly. Also, orbits do not come as close to the horizon for a = 0.8M as they do for a = 0.998M . 
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FIG. 7. Inspiral trajectories in the strong field of a Kerr black hole with a = 0.05M. The top panel includes the effects 

of down-horizon flux; the bottom panel does not. In this case, inspiral is quicker when the horizon flux is included; the hole's 
event horizon acts as an energy sink, as simple intuition suggests it should. This is because the tidal bulge raised on the hole 
by the orbiting body tends to lag, rather than lead, the orbit. Thus the bulge's torque on the orbit opposes the orbital motion, 
causing it to spiral in more quickly. 
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FIG. 8. The + polarization of the gravitational waveform for the inspiral trajectory that begins at t = 40° about the 
a — 0.998Af black hole, viewed in the hole's equatorial plane. The upper panel is the waveform at very early times; the lower 
panel shows the waveform shortly before the inspiraling body plunges into the hole. Notice the very different time scales in 
the upper and lower panels. This is because of the "chirping" evolution of the frequencies fi,^ and Qe — at late times they are 
quite a bit larger than they are early on. Many more orbits per unit time are executed late in the inspiral than early. 
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FIG. 9. The phase ^zmfc and ampMtude Himk for waveform harmonics I = 2, m = 2, and k from to 4, and for spiral into 
a hole with a = 0.998M. These amplitudes correspond to measurement in the hole's equatorial plane. Note that harmonics 
other than A; = become fairly strong at the end of inspiral. The odd behavior in the amplitude for A: = 2 — 4 is because of 
poor computational resolution: the crude parameter space resolution used here was not good enough to accurately capture the 
change in Z^/^ at high frequencies as the body spirals in. 
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FIG. 10. The phase 'i'lmk and amplitude Tiimfc seen in the equatorial plane, for waveform harmonics I = 2, m — 0, and k 
from to 4, and for spiral into a hole with a — 0.998M. In this case, no frequencies are high enough to cause severe resolution 
problems, as was the case for I = 2, m — 2. Note that the phase evolves from high to low late in the inspiral. This is because 
the phase is dominated by the behavior of fie, which chirps backwards when the body gets close to the event horizon; cf. Fig. 

This is essentially a gravitational redshifting effect. Although this is in principle an interesting signature of the black hole's 
strong field, the wave amplitude gets weak as the frequency gets small. This is not surprising, since the rate of change of the 
binary's multipole moments decreases with the frequency. 
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FIG. 11. Evolution of the orbital frequencies fi^ (solid line) and Qe (dotted line) for inspiral beginning at t = 20°. The 
upper panel is for inspiral into a black hole with a — 0.998Af , the lower for inspiral into a hole with a = 0.3594M. In both 
cases, the behavior of fi,/, is qualitatively "normal": it monotonically chirps upward in a very familiar manner. In the CEise 
a = 0.998M, the behavior of Cl$ is comparatively unusual. This is because these orbits are very near to the event horizon 
— their motion in 9 is slowed as seen by distant observers. When o = 0.3594M, no orbits come close enough to the horizon 
for this slowing to have a significant effect. In fact, all orbits in the chosen observation band are at relatively large radius for 
a = 0.3594M, which is why all frequencies are quite smaller than when o = 0.998M (note the different vertical scales in the 
two panels). 
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